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Abstract 


An analysis of the various mechanisms of electromaqnetic wave generation 
by the Shuttie-borne orbiting tether of the T.S.S. Facility has shown that 
significant electrodynamic power levels are available even when overestimating 
the loss mechanisms expected to intervene. For instance, with the tether de- 
ployed downwards, electrodynamic power levels of a few kilowatts are generated, 
when the radius of the terminating balloon is 5 meters or larger. For tetiiers 
deployed upwards (and when an electron gun onboard the Shuttle keeps it at 
plasma potential) electrodynamic power levels of 10 kilowatts are available, 
for balloon radii as small as Z meters. 

This electrodynamic power is in part dissipated by Joule losses in the 
tether. In part goes to accelerate electrons through the sheath surrounding 
the balloon (when in a downward deployment), and in part goes into e.m. wave 
generation. The Alfven Wings are the most interesting wave generation 
mechanism that is activated by the orbiting tether, although from the stand- 
point of practical applications, the use of the tether as a driven antenna, by 
pulsing its DC current, appears the most promising one. A preliminary estimate 
shows that a 100 km tether in orbit would produce ULF/ELF signals that are 
detectable on the ground with state-of-the-art magnetometric instrumentation. 

The power that is expected to go into the acceleration of secondary 
electrons produces effects that range from the excitation of e.m. instabilities 
and consequent e.m. emissions at the gyrofrequency (emissions that are reputed 
to be at a low power level), to the excitation of luminous phenomena, extending 
into the UV band by the bombardment of upper atmospheric layers by the accel- 
erated electrons. 

High-priority topics that require further study are the point above of 
the excitation of instabilities and generation of luminous emissions by accel- 
erated electrons, as well as the generation of ULF/ELF waves by using the tether 
as a driven antenna. 
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1 , Introduction 


1.1 Rationale of the Investigation 

The theoretical investigation illustrated in this report concerns the 
generation and injection of electromagnetic waves in space plasma by means of 
a long-orbiting tether moving in the ionosphere. 

This is a general subject, since there are several electrodynamic 
phenomena associated with the moving tether and, correspondingly, a number of 
problems to be tackled. Before entering a quantitative study for any of these 
phenomena, it is necessary first to have a physical understanding of them, and 
second, to be able to estimate their relative importance. These are in fact 
the basic purposes of this report. The study concludes with a discussion of 
which phenomena are worth deeper study in view of an observational program 
using the Shuttle-based TSS facility. 

1 . 2 The Various Mechanisms of Electrodynamic Interaction and Generation of 

Electromagnetic Waves ' ' 

The configuration that we will study involves a metallic tether connected 
to the Shu .tie at one end, and to a conducting balloon at the other. The tether 
can be deployed both upward and downward with respect to the Shuttle, which 
amounts, from the electrodynamic point of view, to exchanging the polarity of 
two end electrodes (the balloon and the conducting part of the Shuttle) with 
respect to the plasma, as indicated in Figure 1. 

The results we will report refer primarily to a tether moving perpendic- 
ularly to magnetic field lines, at the Shuttle velocity: 

V = 7.8 kra/sec 

Length and radius of the tether are taken as: 

L = 100 km, r^ = 0.5 min 
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BALLOON 



Figure 1. Shuttle-based electrodynamic 
tether 
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A long-conducting tether moving in the earth's ionosphere can lead to 
generation of electromagnetic waves in different frequency bands from ULF to 
MF, through several mechanisms (Williamson and Banks 1976; Dobrowonly et al .. 
1979; Dobrowolny 1979; Banks et al . . 1980), 

(A) Spontaneous Wave Generation 

As the tether system moves through the earth's magnetic field, a 
polarization electric field: 

E = - V X B 

is seen along the tether from the plasma rest frame (we refer, for the moment, 
to the case of a perfectly conducting tether), With a tether length L = 100 km, 
a Shuttle velocity v = 7.8 km/sec, and Earth's magnetic field B = 0.3 gauss, we 
obtain a maximum potential difference between the ends of the system: 

Vq = 2.34 X 10'* volts 

(for y^J. B^). A corresponding potential difference is therefore seen to be 
applied between the lines of forces that pass through the terminations of the 
system and lead to propagation of waves away from the region of the disturbance. 

We can speak of "spontaneous" wave generation in the sense that the 
system is completely passive (no onboard transmitter and no pulsing of its 
natural current) and the original source of energy is provided by the tether's 
motion. As has been pointed out (Gross i and Colombo 1978; Dobrnwolny 1979), 
this electrodynamic interaction (and corresponding wave generation) is similar 
to the interaction of certain celestial bodies with plasmas, like the moon lo 
of Jupiter, in its motion in the Jovian magnetosphere. Recent Voyager I measure- 
ments in the lo's flux tube (Ness et al . , 1979) have actually confirmed the 
occurrence of this interaction. 


3 . 


(B) Driven Wave Generation 

Through proper modulation of an electron emitter at the Orbiter, one 
can obtain pulsating currents in the conducting tether (Banks et al ., 1980), 

which then becomes a long driven antenna. One can also use a transmitter to 

obtain the pulsating tether current. A transmitter wo- 5 ''w necessary in the 
case of a horizontal tether oriented in the direction or inuw there is 

no natural current in the system. 

There would be no dynamical problems with this configuration. In fact, the 

horizontal tether, if terminated at the free end by a balloon (that operates as 

an aerodynamic brake) would be a stable system configuration (Colombo, 1980). 

(C) Generation of Accelerated Electron Beams 

In configuration A of Figure 1 (balloon downward with respect to the 
Shuttle), and under the assumption that the Shuttle is kept at low potentials 
with respect to the medium, large potential drops may be present between the 
balloon surface and the ionospheric plasma. Then, as first pointed out by 
Dobrowolny (1979), secondary electrons can be produced by the impact of 
energetic ions on the surface of the balloon. These secondaries may then be 
accelerated away from the balloon along magnetic lines towards Earth's atmo- 
sphere. [Emission of photoelectrons from the balloon surface is not of primary 
significance at the altitudes under consideration (Dobrowolny 1979)]. Thus, 
generation of accelerated electron beams is another phenomenon that might be 
associated with the tether's interaction with the surrounding medium. This has 
been suggested by a number of authors (Goldreich and Lynden-Bell 1969; Gurnett 
1972; Hubbard et al . , 1976) in studies of the interaction of lo with Jupiter's 
magnetosphere. 

(D) Generation of HF (High Frequency) Electromagnetic Waves 

The generation of accelerated electron beams may give rise to instabili- 
ties, through electron cyclotron resonance interactions, and may produce 
whistlers (or even higher frequency waves),, At the altitudes of interest, the 

electron cyclotron frequency is MHz and the electron plasma frequency 
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fpe ranges from 5 to 11 MHz. It is therefore in the HF (high frequency) band 
where we may expect wave generation by accelerated electron beams. 

1 . 3 Plan of the Investigation 

As the previous summary indicates, there is a rich phenomenology asso- 
ciated with the electrodynar.de tether and, correspondingly, many different 
associated problems fall under the general heading of wave generation. 

Prior to any detailed investigation of the generation of some specific 
waves, it is essential to be able to calculate, for a given tether configura- 
tion, how the primary electrodynamic power available (associated with the 
V X ^ force) is directed into different channels. In particular, we need to 
determine how much power is dissipated into ohmic losses in the tether, how 
much (if any) goes into acceleration of electron beams, and how much is pre- 
sumably left for wave generation. 

We need a model to compute the current-voltage characteristics of the 
long tether (with end electrodes), taking into account finite resistance of 
the tether and, where present, the effect of photoelectric or secondary 
emission. This will be developed in Section 2 of this report. Results will 
be derived for a passive tether (without ion or electron guns on the Shuttle), 
for different tether resistances, and as a function of the radius of the con- 
ducting balloon that terminates the tether at one end. Some comments and 
results on the use of electron and ion guns will also be given. 

In Section 3 we will consider possible effects connected with secondary 
accelerated electrons. A discussion of the distribution of emitted secondaries 
will indicate that the excitation of electromagnetic instabilities is probably 
not important from the energetic point of view. Excitation of UV radiation by 
the high-energy electrons at the foot of the tube cut through by the balloon 
can be an interesting possibility, which is inferred from an analogy with a 
similar phenomenon measured at the foot of tiie flux tube of Jupiter's satellite 
lo. 
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In Section 4 we will discuss the mechanisms of spontaneous wave genera- 
tion and the range of frequencies of the electromagnetic disturbance associated 
with TSS, We will then focus our attention on the generation of low-frequency 
Alfven waves, which requires large balloon dimensions. The physical picture of 
the Alfven wave system associated v/ith the TSS motion will be described. 

Current along magnetic flux tubes, power, and impedance associated with trans- 
mission of Alfven waves from the tether will be computed. The Alfvenic current, 
obtained as a function of balloon radius, will then be compared with the current 
obtained from the model for particles collected from the plasma by the end 
electrodes. 

In Section 6, on the basis of the calculation of power in Alfven waves, 
we will use theory of wave propagation to the ground for a point dipole source 
in the ionosphere (despite its inadequacy for our system) to determine the 
potential for detecting VLF waves at ground level, and will draw a preliminary 
positive conclusion. Finally, Section 5 will summarize the main points of 
the report. 


2. Studies of Current-Voltage Characteristics of the Flectrodvnamic Tether 
2.1 Model for Charged Particle Collection 

The tether-balloon system, moving across magnetic lines, carries a current 
due to the collection of ionospheric ions and electrons at its ends (the balloon 
and the conducting part of the Shuttle, which is assumed to have an area of 35 m^). 

The condition that one imposes to derive current and potentials is that of 
a balance of charged particle fluxes between the two end electrodes. In terms 
of currents, the current collected at the Shuttle i^ has to be equal and 
opposite to the current ig collected at the balloon: 

= -ieUg) (2.1) 
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(for a tether system without any gun, ion or electron, at the Shuttle). In 
(2.1) we have explicitly indicated that the currents are functions of the 
potentials (Vg,Vg) of the two electrodes, with respect to the plasma. For a 
perfectly conducting tether, it would bes 

IVs-Vjl - |Vo X B • L| (2.2) 


where v., 1s the Shuttle velocity (v„ = 7.8 kin/sec), B Is Earth's magnetic 
field (B 0.3 gauss, at the altitudes of interest between 100 and 300 km), 
and L is the tether's length. Hence (2.1) is an implicit equation for one of 
the potentials, for example Vg. Having determined Vg,Vg is obtained from (2.2), 
and then the current 1^. = |ig| = |igl is calculated from the solutions found for 
the potentials. 

The model that has been used for charged particle collection is the follow- 
ing (Anderson e t al . . 1979). For the attracted particle contribution to the 
current referring to particles of species j (j = i, e for ions and electrons 
respectively), we write: 


i. attracted 

w 


i 


Jo 


= f (V) 


(2.3) 


where the normalizing electron and ion currents are given by: 

''lo “ 1 V 


(2.4) 

(2.5) 


(with A the collecting area, the electron density, and v^^,g thermal velocity). 
The function f, which depends from the electrode potential V, through*. 


'kTp 


(^) 


4/3 


(2.6) 


(Tg being electron temperature, X^g the Debye length, and R the electrode's 
radius), is plotted in Figure 2. The function f was derived by combining di f- 
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ferent models for particle attraction by large electrodes at large and moderate 
potentials (Alpert et al . . 1965; Linson 1969). 

For the repelled particle contribution to the current, we used simply; 



where V Is now the repelling potential. 

Approximate results for currents and potentials as a function of balloon 
radius, obtained with the above model for a perfectly conducting tether, are 
given by Anderson et al . (1979). 

Recently, we have developed (1980) a more accurate method, based on a 
transmission line analogy of the tether system, to compute the stationary state 
described by (2.1). We did this by solving a time dependent problem, and hence 
obtained also the transient of the tether system towards the stationary state. 

This method, which was devised especially for the more difficult case of 
computing current and potential distributions of a bare metallic tether, has 
now been applied to the tether-balloon system under consideration (conducting 
insulated tether with terminal electrodes), including also the effect of tether 
resistance and, as will be described, the possible emission of secondary 
electrons from the lower electrode of the system, and the effects of using 
electron or ion emitters. 

Finally, it should be noticed that the equilibrium current i obtained in 
the tether through the above mentioned methods is obviously a very non linear 
function of the equilibrium potential difference (Vg-Vg), also obtained from the 
calculations. The various results for current and potentials versus balloon 
radius which will be seen in what follows, do not correspond therefore to a 
ohmic relation between i and V^-Vg involving the tether resistance It is 
however tre that, in all cases, the currents obtained must be smaller than the 
ohmic limit that corresponds to the total electromotive potential available, i .e . : 
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2.2 Results for Current dnd Voltages 

We first present results for current 1n the tether and voltages of the 
Shuttle and the balloon for a completely passive system— I.e .. a system not 
using Ion or electron emitters, Figure 3 gives the variations of the Shuttle 
(Vg) and balloon (Vg) potentials as a function of balloon radius Rg for con- 
figuration A (tether deployed downwards) and a resistivity of the tether p ■ 
0,03ygm (aluminum- type tethers). In Figure 4 we have the same quantities for 
a case of higher resistivity, p * O.lBpom (piano-steel -wire tether). As can 
be seen essentially In both cases, the Shuttle potential Vg, which Increases 
with balloon radius, reaches a/IO^ volts at Rg 'v 10m, At higher balloon radii, 
Figure 4, which refers to higher resistance, gives somewhat smaller values for 
Vg than for the low resistance case of Figure 3, 

Another qualitative difference between the results of Figure 3 and 
Figure 4 is given by the drastic drop in balloon potential |Vg| which occurs 
at Rg = 85 - 90 m for p « Q.lSpC^i (see Figure 4), whereas for the lower 
resistivity p = 0.03pom, lVg| decreases with increasing radius but remains 
relatively high up to the maximum radius (Rg » 100 m) that we have calculated. 
Clearly, the Shuttle potential must increase with radius in order to balance 
the increasing current collected by the balloon; eventually, the balloon 
potential must drop to very low (negative) values. The reason why it does so 
at smaller radii for the case of the Higher resistance depends only upon the 
increased ohmic losses (at any given radius) that correspond to such higher 
resistances, 

At large balloon radi<, as could be necessary to generate larger cur- 
rents and reach substantial power levels in Alfven waves guided by the balloon 
flux tube (see Section 4), the potential of the Shuttle in such a passive 
system can be relatively high. The use of an ion emitter at the Orbiter would, 
therefore, be necessary to keep its potential at low levels. 

Figure 5 gives the variation of current in the tether as a function of 
balloon radius, again for configuration A and the two values of resistivity 
considered before. 



The difference between the two curves of Figure 5 and, In particular, 

the dip of curve b) (p “ O.lBunm) at Rg * 85 ni can be explained In terms of 

the different behavior of the potential jVgl as a function of the radius for 

the two cases considered (see Figures 3 and 4). In the case of the higher 

resistance, at a radius between 80 and 90 ni, the current reaches a quasl- 

plateau. We understand that on the basis of the formulas for collection 

efficiency given In Section 2.1, the balloon current Is proportional to Rg 

f (|Vgj). In the region (see Figure 4) where |Vgl drops to small values, 

the collection efficiency f drops correspondingly. This roughly balances the 

Increase of current with Rg and produces the plateau. Once lVg| has dropped 

to small values (l|-j^| i 1)> f (|Vgl) is essentially 1 (see Figure 2) and 

^ 2 

the current starts incru'slng again (as Rg), This explains the dip of curve 

b) In Figure 5. The same effect Is not present in curve a) because, for the 

case of the smaller resistivity (p = O.OSpnm) there Is no drop In |Vg| for 
Rg s 100 m, as seen in Figure 3. 

The same quantities ( 1 . e . , voltages and current) are plotted In Figures 6 
and 7 respectively, referring now to configuration B (tether deployed upward) 
and for a resistivity p = O.lSynm. The curves for the smaller resistance are 
not very ,'lifferent. These results are interesting because they show that, 
even with small radii of the balloon (or subsatellite)-- e.g., Rg 'v Im— in a 
passive system, the Shuttle goes to high potentials with respect to -the plamsa 
(|Vgl 10^ volts), while the corresponding current Is quite small (1 'v 0.08 
amps). The physical reason for this is the fact that the balloon collects now 
electrons. In order to balance this electron current, the Shuttle, which 
collects the much less mobile ions with its collecting area of 35 m^, must 
necessarily charge to a high potential. Clearly, the use of an electron gun 
at the Shuttle is strictly necessary for configuration B-- both to keep the 
Shuttle potential low and to obtain larger currents. 

With either configuration (and also with or without the effect of 

secondary electron emission), the primary limitation on tether current is seen 

to be the resistance of the tether. For the two values of resistivity 

12 , 
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considered, the limiting values are: 

i|^ = 1.23 amps for p “ O.lBpfim 
ij^ = 6.1 amps for p = O.OSpQm 

From the curves of Figures 5 and 7, we see that these limiting values are not 
reached for a passive system, even at large balloon radii. Such values will, 
therefore, be approached only through the use of charged particle guns. 

2 . 3 Electrical Currents with the use of Charged Particle Guns at the Orbiter 

Referring to configuration A, suppose we want to use an ion gun at the 
Shuttle in such a way as to keep its potential with respect to the plasma at 
very small values (Vg'^O). Then the balloon will be, approximately, at a 
potential - jVgj where: 

IVgl = - R^i (2.8) 

and - |VqXB_*L|. current i will be given by (see Section 2.1): 

= ■'b = 'io ^ (2.9) 

For a given balloon radius Rg, the actual values of current in the tether and 
potential of the balloon are obtained from the intersection of the two curves 
i(Vg) given by (2.0) and (2.9). 

The result for the current as a function of balloon radius is given in 
Figure 8. The Interpretation is that, for a given balloon radius Rg, an ion 
gun at the Shuttle with current i(Rg) obtained from Figure 8 is necessary in 
order to keep the Shuttle at essentially zero potential. The curve, which 
refers to the low resistivity tether (p = 0.03ynm), gives obviously larger 
currents than those obtained for configuration A and a purely passive system 
(see Figure 5). 

In the same- way we have obtained in Figure 9 the current of an electron 
gun necessary, in configuration B, to keep the Shuttle potential at low levels. 
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Comparing with Figure 7^ referring to a passive system, we see that the use 
of an electron gun for configuration B, besides being necessary to keep the 
Shuttle potential low, is the only way to obtain significant current values 
regarding electromagnetic wave generation (for example, using the long tether 
as a dipole antenna). For balloon radii Rg>10m, one already approaches the 
resistive limitation (ij^'vB.l amps). 

2 . 4 Emission of Secondary Electrons by Pr oto n Impact on the Balloon Surface 

As can be seen from Figures 3 and 4, referring to a passive system, the 
downward balloon (configuration A) can be charged to a very high and negative 
potential. This is, of course, true also if we use an ion gun at the Shuttle 
to keep its potential low, unless the current of the gun approaches the resis- 
tive limit. 

In these conditions, the ions impacting on the balloon surface can reach 
high energies through the accelerating sheath region surrounding the balloon, 
and then can cause significant emission of secondary electrons. Since this 
effect will vary the current, and hence the voltages, in the system, it has to 
be taken into account in computing current-voltage characteristics. The effect 
is also of interest because the secondary electrons, accelerated away from the 
balloon through the sheath region, could produce electromagnetic waves (as 
mentioned in the Introduction). 

We will come back to this second aspect later. For now we will focus 
on how to take the effect of secondary emission into account in our calculation 
of current-voltage characteristics. 

For a given material of the balloon surface, the yield of secondary 
electrons per ion impact will depend on the energy of the impacting ion and 
on its angle of impact, with respect to the normal to the surface (Massey and 
Burhop 1952, McDaniel 1964). In order to discuss how this energy depends 
upon the potential Vg of the balloon with respect to the plasma, we have to 
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consider the ion motion in the sheath region. First of all, we can suppose 
the sheath electric field to be radial; this is justified by the fact that wave 
effects are probably not important as: 


1 

2 



« leVg 


(im|v^ 'V' 5eV, jeVgl -v 2.3 x lO** eV). 

For high balloon potentials, we can write for the sheath thickness 
(Alpert 1965, Dobrowolny 1979): 


*^sheath 

Rb 


<v 0,8 



( 2 . 10 ) 


Using parameter values for an altitude of approximately h = 120 km (the 
altitude of the balloon) and using the results for the potential |Vgl shown 
in Figure 3 (for p = O.OSpnm), we obtain for the thickness of the sheath, as a 
function of balloon radius, the values reported in Table I. For comparison, 
the ion Larmor radius (0^ ions) and the electron Larmor radius are given by: 

R^- -V 4.7 m , Rg 2.7 X 10"^ m 


If we now go to the equation of motion: 

mi ^ = q [E (r) + V X B] 


(2.11) 


and compare the electrical and Lorentz force, we have, in order of magnitude: 


where: 



a 


1 2 
— m.v 
? 1 max 



( 2 . 12 ) 


(2.13) 
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Is the maximum energy the ions can pick up through the sheath region. 
Estimating? 


E -v ^ 
Rsheath 


we obtain: 




^*^sheath 


VJ/2 
-V 9.6 J. 

Rsheath 


(2.14) 


The values of a corresponding to various values of balloon radii are also re- 
ported in Table I. It is a > 1 and, therefore, the electrical force dominates 
in the equation of motion (2,11). In these conditions, we have a radially 
accelerated motion and all the ions entering the sheath region will finally 
pick up the energy (2,13) corresponding to the potential drop [Vgl between the 
balloon and the plasma (and all of them will normally hit the surface). 

Hence, in the stationary state, the current balance condition (2.1) (for 
configuration A) can now be rewritten, taking secondary emission into account 
as; 

i = (V^) = (V^) [1 + S (|Vg|)] (2.15) 

where the yield 6 is a function of the total balloon potential with respect to 
the plasma. 


Table I 

Values of the Parameter a 



1 I 

*^sheath 

a 

(meters) 

(volts) 

(meters) 

1 

2.33 X 10"^ 

7.1 

2.07 X 10^ 

5 

2.29 X lO'* 

14,07 

1.04 X 102 

10 

2.24 X lO** 

18.74 

7.69 X lOi 

20 

• 2.11 X 10"* 

24.59 

5.69 X 10^ 

30 

1.96 X 10*+ 

28.34 

4.76 X 101 

50 

1.69 X lO't 

33.15 

3.78 X lOi 
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2.5 Discussion of a Model for the Yield in Secondary Electrons 


The main problem now Is the choice of the function 5(Vg), To our 
knowledge, there are essentially no data on secondary emission from Impacting 
o’** Ions (the dominant species In the Ionosphere at the altitudes of Interest). 
Besides, the data available are, In general, very scarce for any 1on (not only 
o’*”) In the range of energies from 1 to 20 keV, where the effect of secondary 
emission could be significant. 

Figure 10 gives the yield as a function of energy for protons impact- 
ing on two different materials. The yield was obtained by Whipple (1965) by 
putting together data from different experiments, at various energies. Includ- 
ing the range of Interest here. Although we do not have an analogous curve 
for 0 ions, we can still derive from the existing literature on secondary 
phenomena (McDaniel 1964) some useful Information on the variation of the 
yield, with Impact on mass at various energies. 

What one finds Is that, at low energies (<1 keV, approximately), the 
yield decreases with mass (see curves referring to ions of various noble 
gases by Hagstrum 1956), On the other hand, for very high energies (>100 keV), 
the yield seems to be almost proportional to the mass (Massey and Burhop 1952, 
McDaniel 1964), There are almost no results for Intermediate keV energies, 
where the two dependencies from mass (that appropriate for very low energies 
and that for very high energies) should obviously merge into one another. 

Figure 11, taken from Higetsherger et al . (1956), goes up to 8 keV, includes 
o’*’ Ions, and Indicates a smaller yield for Increasing mass. 

On the basis of this Information', we will use In the following study of 
current-voltage characteristics the Whipple's curve of Figure 10, referring to 
protons (and the aluminum target). The discussion Indicates that, in the range 
1 keV - 20 keV, this curve may give an overestimate. For energies lower than 
'v/l keV, where the yield decreases strongly with mass, the use of the proton 
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Figure 10. Secondary yield of electrons for 
proton impact at energies above 
L KeV (from Whipple, 1965). 
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Figure 11. Electron yields of various Ions on 

Ag-Mg targets (from M.J. Hlgatsbergerj. 
H.L, Demorest, and A.O. Nier, J. ApdI. 
Phys . 25 , 883, 1854), 
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curve certainly re, /esents an overestimate, un the other nano, ror these 
potential values, the yield (for protons as well) Is already quite small and 
does not have any appreciable effect on current-voltage characteristics. 

2.6 Results for Current and Voltages In Presence of Secondary Electron Emission 

Figures 12 and 13 give the Shuttle and balloon voltages as a function 
of balloon radius, with secondary electron emission taken Into account. They 
refer respectively to the two values p » O.OSuam and p - 0.15u«m of the 
resistivity of the tether. They should be compared with the curves of Figures 
3 and 4, referring also to a passive system* and configuration A, but not taking 
the effect of secondary emission Into account. 

Figure 14 gives current as a function of balloon radius for the two 
previous values of resistivity, and includes the effect of secondary electrons. 

Emission of secondaries increases the current, as can be seen by comparing 
the curves of Figure 14 with the corresponding curves of Figure 5, For example, 
for p « 0.03pfim and Rg « W m, the current goes from 0.23 amps to 0,98 amps 
when secondary emission is taken into account. The increased current, which 
leads to greater ohmic losses, tends to lower the balloon potential and hence 
transfer some of the potential drop to the Shuttle. This effect cannot, how- 
ever, be very large; in such a case, the emission of secondaries, and hence 
the current, would start decreasing. Comparing potential values for correspond- 
ing balloon radii between Figures 3-4 and 12-13, we can see that there is in 
fact some displacement of potentials in the sense mentioned above. In spite 
of the fact that the effect of secondaries is probably overestimated, these 
variations are nonetheless rather small. 

Notice that the curve b) in Figure 14 (that refers to p = O.lBpnm) 
shows a monotonic behavior contrary to the corresponding curve in Figure 5. 

The reason for the difference is that, when the balloon potential drops to very 
low values (that explain the feature of curve b) in Figure 5), the contribution 
to current due to secondary electrons drops to zero. The change in slope of 

curve b) in Figure 5 does not appear in Figure 14 because of the effect above. 
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Figure 12. Configuration A : Voltages versus R 
taking secondary electron emission 
into account for L= 10^ m, =0.03 uilm. 
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Figure 13. Voltages versus R„, taking into account 

Ij’ 

secondary electron emission, Configurat] 
L = 10^ m, M = 0.15 jdSlm. 


27 . 




Rp (meters) 


Figui-G 14. Configuiratloii A : current versus Rg for L » 
and taking secondary electron emissions into 
account. 
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2.7 Partition of Primary Electrodynamlc Power into Different Loads 


Prior to detailed investigation of any specific mechanisms or effect of 
the tether's electrodynamlc interactions, it is of fundamental importance to 
understand how the primary electrodynamic power of the tether system (coming 
from the y. X B electric field) is divided into different loads. 

The power that we calculate from: 

Pi = IVg-Vji (2.16) 

by using the results obtained for current and potentials, does contain the 
effect of ohmic losses in the tether (i .e. . |Vg-V[^| = V^-R^^i). The ohmically 
dissipated power is given by: 

Pz = R„1" (2.17) 

If there is no effect of secondary electron emission, we can say that the 
total, power P 4 available for exciting wave processes outside the tether is 
equal to Pi. That is: 

P 4 = Pi (2.18) 

Where acceleration of secondary electrons takes place, an overestimate for 
the power P 3 in accelerated electrons is obtained by supposing that, in pass- 
ing through the sheath surrounding the balloon, the electrons gain all the 
energy |eVg|. Then: 

P 3 = IVjl 133, (2.19) 

with igg^ the current in secondary electrons. This is indeed an overestimate, 
partly because the emitted secondaries are accelerated through the sheath 
region only through the component of the electric field parallel to which 
adds a numerical factor (<1) to the above formula. Since our treatment of 
the secondary emission is an approximate one (see Section ,2.5), and since the 
current is also overestimated, we will use (2,19) as it is. There: 
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sec 



(Vb) (iVBl) 


At tliQ oqu1Hbr1uni state we can use (2.15) to obtain, in terms of the total 
current i: 


i 


sec 


« i 


■S(Vb) 

1 + 6 (Vp) 


( 2 . 21 ) 


and, correspondingly: 


Pa 


s 

1+6 


( 2 . 22 ) 


In the general case P 3 0, the power Pi,, which is available for wave genera 
tion, is given by: 


Pi, ° Pi - P3 


(2.23) 


More precisely, Pi, is power going into wave generation (once, for example, 
the tether current is pulsated), if there are no other dissipative phenomena. 
Generation of microturbulence in the charged sheath surrounding the balloon 
could absorb power Thus, Pi, must be considered to be rather an upper limit 
for the power available to wave generation. 

2 . 8 Results for Different Power Levels 

Figures 15 through 18 give the results for the computation of the power 
levels P 2 , P 3 , and p 4 as a function of balloon radius for a passive system 
(configuration A) and both with and without the acceleration of secondary 
electrons. Again, the two values of resistivity (p - O.OSynm and p = O.lSynm) 
are considered. 

The following important points can be made from these results: 

(a) Even with high values of current (obtained for large values of Rg) 
and in spite of the high ohmic losses, the power Pt, available for wave genera- 
tion is substantial, of the order of several kilowatts (even for the case of 
higher resistance; see Figures 17 and 18). 
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Figure 15. Power into ohmic dissipation (P 2 ) and into wave 
radiation (P 4 ) without effects of secondary 
electrons. Configuration A, L = lO^m, ^=0.03 pilt 
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Figure 17. Power Levels versus R 
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(b) The power P 3 , which goes into acceleration of secondary electrons 
is also significant. It amounts to several kilowatts— even at very large 
radii for the low resistance case (see Figure 16)— whereas it drops to lower 
values for larger radii, due to the higher ohmic losses, in the case of higher 
resistance (see Figure 18). Even in this last case, P 3 > 10^ watts up to 
radii Rg < 30 m. Although this power, as it comes from the previous discussion, 
is certainly overestimated, its values make it potentially interesting to dis- 
cuss effects produced in the medium by the accelerated secondary electrons. 

(c) If we compare the wave power P 4 , obtained with and without secondary 
emission effects, we see that the effects are not much different. 

Substantial power levels are therefore available for wave generation, 
even if the effects of secondary electron production and acceleration are 
taken into account. This important conclusion would stand even in a more exact 
treatment of the problem, because it comes from an overestimate of the effect 
of secondary electrons. It is clearly a very significant conclusion in terms 
of the capability of the long tether to radiate electromagnetic waves. 

Besides these general conclusions, it is worth to comment on the be- 
havior of the curves P 4 that represent the power available for wave radiation, 
in Figure 17 and Figure 18 (p = O.lSynm). In these Figures, we have performed 
calculations using larger balloon radii than in the previous Figures. In both 
cases there is a decrease in Pi* starting with a radius of about 90 m, while P*,, 
and the ohmic power loss P^ as well, go to an asymptotic value for radii 
Rg s 100 m. This can be understood in terms of the current balance that we 
are imposing, and of the formulas for current collection that were given in 
Section 2.1. A consequence of these is that both IVg-Vgj and the current i 
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go necessarily to asymptotic values for large enough values of Rg. What 
happens is that, for p ** O.lByjmi and Rg » 90 m, the balloon potential |Vg| 
drops to very low and negative values (see Figure 4), If jVg| continues to 
drop from these low values, some electrons start to be collected at the balloon 
(in spite of its repulsion potential), This causes a decrease in current and 
explains the decrease in Pi, around the value Rg » 95 m, as seen in Figure 17 
and 18. However, if |Vgl becomes too small, the current eventually decreases 

so much that the ohmic losses decrease and hence |Vg-Vgl starts increasing. 

This means that |Vg| increases slightly again and settles both IVg-Vgl and i 
to their final asymptotic values. Hence, the same happens for Pi, and ?z- 

In order to observe the same behavior on the curves of Figure 15 and 
Figure 16, that relate to the case of lower resistivity, we would have to go 
to much higher values of Rg, because this effect is critically dependent upon 
ohmic losses, which are now smaller. 

Calculations of powers referring to configuration A (with an ion gun at 
the Shuttle such that V^'^O) will not be reported here, since they leave the 
previous general discussion unchanged. 

We report instead, in Figure 19, the power P 4 available for wave genera- 
tion, versus balloon radius for the more practical case of configuration B 
(balloon upward) with an electron gun at the Shuttle keepi“'g the Shuttle' 
potential low. (The calculation refers to p = O.OSyfim.) The values of power 
are substantial: above 10 kilowatts up to R0'^3 m. They then decrease with 
radius on account of the high currents (see Figure 9) and the correspondingly 
high ohmic losses, but still remain substantial for the radii considered. 
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3. Effects of Accelerated Secondanv Electrons 


3.1 Introduction 

The Investigation of possible processes generated by accelerated secondary 
electrons poses a difficult problem. It requires, first of all, an Investigation 
of electron trajectories In the sheath region surrounding the balloon, to 
determine their distributive function outside this region, Second, one has to 
Investigate the stability of such distributions in the background plasma with 
respect to the generation of electromagnetic waves, and actually go to the non- 
linear stage of Instability In order to be able to calculate the efficiency of 
these processes. A final problem concerns those high-energy electrons that have 
possibly survived any of these interactions: we must evaluate their energy 
degradation in the atmosphere and the possible resulting excitation of luminous 
phenomena. We will not go Into all these detailed mechanisms here, Consistent 
with the rest of the report, however, we will outline the various phenomena 
involved to get an Idea of their possible relevance for future studies. 

3.2 Distribution of Secondary Electrons 

As we have already metioned (see Section 2.4), one may suppose that 
the sheath electric field around the balloon is radial and, in order of 
magnitude, given by: 

^ ''B^''sheath 

Values of the sheath thickness have been given for various balloon radii in 
Table 1 of Section 2.4. If we consider the electron equation of motion in the 
sheath region: 

mg = e [E + V X B] (3.1) 
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and comparG electrical and magnetic force (as we have done for the ions In 
Section 2,4), we are led to define the parameter: 


a 


(e) 



(3.2) 


wi th ; 


1 

2 


111 , 


max 


eV, 


It i 


■>. 5.62 X TO-2 


Vb V2 


'sh 


(3.3) 


and we have reported values of this parameter, for various radii Rg, in Table 
II, having used for Vg the results of Figure 3. Contrary to the ion case, it 
oig ^ i. This means we cannot assume that the electric force is dominating. 
In Figure 20, which depicts the sheath region surrounding the balloon, we have 
Rp “ Rg + Rg^. The electrons will exit the sheath region at radial distance 
r « R from the balloon center and we want to have some idea of their distri- 
bution in a point on the flux tube at distance p from the line of force passing 
through the center of the balloon, 


Table II 

Values of the Parameter 


(meters) 

(volts) 

^sheath 

(meters) 

a(e) 

1 

2,33 X lO** 

7.10 

1.208 

5 

2.29 X 10'* 

14.07 

0.604 

10 

2.24 X 10‘* 

18.74 

0.449 

20 

2.11 X lO't 

24.59 

0.332 

30 . 

1.96 X 10‘^ 

28.34 

0.277 

50 

1.69 X lO'^ 

33.15 

0.220 
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within the sheath, the electrons are accelerated parallel to B by the 
parallel electric field component; 

Ell « cosO (3,4) 


and, with the notation of the figure, 


slnO « (3,S) 

Rr 

The velocity Vn^ (p), which they attain after tranversing the sheath under 
the accelerating force. Is given by*. 


v,,„ (p) = (1 - 

0 


(3.6) 


It tends to 0 forfi..+ 1, such as on the boundary of the flux tube; the electri 

^r 

field (radial) tends to be perpendicular to B. 



Figure 20. Geometry of the sheath region 
surrounding the balloon 
■ 39. 



The motion perpendicular to B is, to a good approximation (as the electron 
Larmor radius Rg is « Rg|^), given by the E. x J| drift and the corresponding 
perpendicular velocity is: 

''lo (p) 

This tends to 0 as 0, v/hile the electric field tends to become parallel 

Kr 

to B, The electron distribution function then can be written, for example, 
as; 


(P,, 


VI v„) « nt 


(()) 


exp 


1 ^- 


Vu--VMo(q) 


'the 


•) - (- 


'^the 


Vie 


(3.8) 


where (p^) is the density of the secondary electrons and v^j^^ ('x- 1.5 x 10^ 
km/sec) their small thermal velocity. 


3. 3 Electromagnetic Instabilities ,of Secondary biactron Beams 

To be excited, all electromagnetic modes require an excess of perpendic- 
ular (over parallel) momentum (Melrose 1973, Mangeney and Veltri 1976). A 
practice' condition for any of the electromagnetic instabilities is that: 

V|„ 

y = arct < 25° (3.9) 

''iio 

If we express this condition using the above formulae (3.6) and (3.7) for 
ViiQ vio» we obtain: 

§*> 0.955 (3.10) 

Kr 

or, equivalently, for the angle e of Figure 2, e > 73°. Therefore, possible 
electromagnetic instabilities are expected only quite close to the boundaries 
of the flux tube cut through by the balloon. 

Next, we wil.l briefly discuss the type of electromagnetic instabilities 
we have to expect. Let's look just at the cyclotron resonance condition; 
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V" „ 

w nwcose — • 'u o 

ce c 

kc 

with n = ■^ the index of refraction, o being now the angle between k and B 
and V = 0, ± 1 , ± 2, etc. For the ionospheric condition: 

'v<10 (3,12) 

“ce 

the high frequency modes (ordinary and extraordinary modes with w > w and 

pe 

n < 1) will be in resonance only for v » 1 and, correspondingly, have very 
reduced growth. The only significant instability of electromagnetic waves 
that would remain is that of the whistler mode at w a , 

3 . 4 Significance of Instabilities and Generation of UV Emissions 

Condition (3.8), limiting the generation of whistler waves to the 
boundary of the flux tube cut through by the balloon, is quite restrictive. 

If we take the total power in accelerated secondary electrons (see Section 2.8) 
and take into account the tiny fraction of the flux tube cross-section that is 
significant for whistler wave generation, it is easy to see that the available 
power is reduced to insignificant levels. We could, tnerefore, conclude that 
electromagnetic wave emissions from secondary electrons are probably not of 
importance. 

On the other hand, if we consider the electrons in the remaining part 
of the flux tube cut through by the balloon and, in particular, close to the 
center (which tend to have an excess in parallel velocity), it is conceivable 
that they might excite an' instability of electrostatic plasma waves. A 
quantitative appreciation of this effect could be a topic for future studies. 

In the same way, it would be worthwhile to investigate the possibility 
that accelerated electrons, in the absence of other interactions, reach the 
earth's atmosphere' and finally lead to the excitation of some luminous 
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phenomena (for example, excitation of UV radiation). An appreciation of 
this possibility comes indirectly from a comparison with the observation 
of Lyman a UV excitation at the foot of the flux tube of Jupiter's satellite 
lo from the Copernicus telescope (Atreya et al . . 1977). These authors 
estimate that the interaction of keV electrons of 100 eras ciif^sec"^ with 
the Jovian upper atmosphere produces approximately 100 kR (kiloroenpten) of 
hydrogen Lyman-ct at each foot of lo's flux tube. 

In our case, froM Figure 18, we have for Rg = 5m: 

P 3 'V. 2.5 x 10^ watts 
and correspondingly, a power density: 

PD a- 3 X 10*' ergs ciir^sec~^ 

Although this is an overestimate (and we also have to take into account other 
possible losses due to instability excitation), the number obtained is above 
the value quoted by Atreya et al . (1977). A comparison with lo suggests, 
therefore, that in the case of downward deployment of the balloon, we can 
expect to see UV excitation at each foot of the balloon's flux tube in the 
earth's ionosphere. This observation will require UV instrumentation in orbit. 

4. Alfven Wings 

^ 1 Discussion of the Electromagnetic Disturbance Associated with TSS 

The concept of Alfven wings was introduced in the earlier study by 
Drell qt , a| . (1965). It has gained recent experimental support from the Voyager 
I observations of magnetic field perturbations associated with the flux tube of 
lo in Jupiter's magnetosphere (Ness et al . , 1979) and from further work that 
followed the analysis of observations (Neubauer li'^O). Figure 21, taken from 
Banks et al . (1980), gives a pictorial view of the current wings associated 
with a TSS moving in Earth's ionosphere. It is important to recall some basic 
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Figure 21. 

Schanatic view of the upper and lower cxarrent sheets 
which spread out fron the electrod'/namic tether svstem. 

Ihe periodic darkened regions represent the outavard 
propagation of ww frequency' ALfven waves along the magnetic 
field. Stiere is a net positive charge excess on the 
top wing and a net negative charge density on the laver 
wing . 

(from Banks et al . , 1980). 


43 . 


features of the disturbance associated with the TSS. 

Since the ionospheric conductivity parallel to the Earth's magnetic 
field is extremely large at altitudes above the E layer (and much larger than 
the transverse conductivity), the magnetic field lines can be regarded as 
equi potentials. The ionospheric state is perturbed by the motion of the 
tether (or any large conductor) across magnetic lines. From the rest frame 
of the plasma one sees a polarization electric field: 

E-^y^xB (4.1) 

if we refer, for the moment, to the case of a perfectly conducting tether. A 
corresponding potential difference is, therefore, seen between the lines of 
force intercepted by the ends of the system. 

This perturbed state tries to readjust itself (to the previous equilibrium 
state, with no potential difference across field lines) through the propagation 
of waves from the region of disturbance (mechanism of "spontaneous" wave 
generation— see Section 1.2). These waves, and associated currents parallel 
to £ lines, are carrying away the applied potential differences, or the 
equivalent transverse space charge. 

The problem of the electromagnetic disturbance associated with the moving 
tether can be formally set up as a problem of radiation from a current source 
( i . e . , the classical problem of antenna theory). In the case of radiation in 
a magnetized plasma, by combining Maxwell's equations and Fourier transforming 
in space and time, we obtain the following equation for the space-time Fourier 
transform of the radiated electric field: 

. (k,o)) E. (k,m) = - M j.o (k,a)) (4.2) 

where w, _k are angular frequency and wave number of the radiated waves, 
represents the current source (current density) and the tensor A. • is defined by: 

A., (kio) = n^ (x.x. - 6. J - e.. (4.3) 

ij “ 1 a TJ u 
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with! 

n > if (4.4) 

la 

|( 

the refraction index, x = ffT ®ij dielectric tensor of 

the magnetized plasma. As Stix calculated (1962): 

A s det A.. = 0 (4.5) 

' J 

gives the wave dispersion relation. 

For an observer at rest with respect to magnetic flux tubes, the tether 
represents a moving current source, so that we can write for the current 
density: 

(x, y - Vgt, z) (4.6) 

(For motion in the y direction, see Figure 22.) Consequently, by Fourier 
transforming in space and time and by substituting into (3.1), we obtain: 

(k.w) Ej (i.«) = - i 6 (« - J? (k) (4.7) 


The purpose of writing this equation is to point out that it constrains the 
angular frequency of the radiated v/aves, which must satisfy: 

« ■ ( 4 . 8 ) 


This does not, of course, fix the frequency that depends upon k^, which, in 
turn, depends upon the function (k.) and upon the role played by the plasma 
dispersion in solution (3.7). 

It is quite natural to estimate: 


k 


y 



(4.9) 


where d^ is the conductor's dimension in the direction of motion. This then 
gives a frequency: 

f* = ^ (4.10) 

y 
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Figure 22 . Geometry of the tether 
balloon system. 
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It is important, however, to realize that -this is not the typical frequency 
of the electromagnetic disturbance, but rather must be interpreted as an upper 
limit to the frequencies contained in such disturbance. This is equivalent 
to saying that, for an observation site on a magnetic flux tube, the distur- 
ance contains all frequecies: 

f = ^<f* (4.11) 

while the power radiated in frequencies f » f* will be negligible. 

Only if; 

f* < f . 

Cl 

as it is the case for lo, does all the power of the disturbance go into hydro- 

s 

magnetic waves and, in particular, into Alfven waves, which are then guided 
along magnetic lines. 

In the case of the TSS, this can occur for the disturbance associated 
with the balloon if its radius Rg is greater than a critical value Rg*: 

Rg 5: Rg* 'V' (20 sina) meters (4.12) 

where a is the angle between velocity and magnetic field. If we refer now 
to the tether's cross-section (d 1mm), the frequency on the other hand, 
is: 

^*tether ^ 

Thus, the overall disturbance associated with the TSS will generally contain 
frequencies up into the Megahertz range. Of these wave components, only the 
low frequency Alfven waves strictly propagate their energy (at a velocity v^) 
within the related magnetic flux tubes. The higher frequency components of 
the disturbance have different dispersion properties (Stix 1962). 
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4.2 Parallel Current Associated with Alfven Waves 


The coupling of the tether's system with the plasma medium and, 
possibly, with the lower layers of the ionosphere, occurs through the radiation 
of waves. In particular, if a curent J„, parallel to magnetic field lines, 
is associated with these waves, it continues the tether current into the iono- 
sphere— down to E layer altitudes where perpendicular current closure can 
take place, 

This parallel wave current is beyond the so-called "dc current model" 
of lo, proposed by Goldreich and Lynden-Bell (1969). According to the model, 
the flux tube Intercepted by lo would actually be frozen to the satellite and 
follow its motion around Jupiter, with up-going and down-going parallel 
currents at the boundaries of the tube (the Alfvenic currents) and circuH 
closure-- within lo on one side, and across Jupiter's ionosphere on the other 
side. 

We will now derive a general equation, with no approximation of small 
amplitude for the waves, for the parallel current associated with Alfven 
waves. The Alfven waves characteristic's are given by (Jeffrey and Taniuti 
1964): 

B 

V ± T = constant (4.14) 

- 

where y. and ^ refer to the fluctuations of velocity and magnetic field in the 
waves and ppis the plasma mass density. The constant can be evaluated from 
the background properties of the plasma, In the rest frame of the plasma, 
then; 


V £ 




1/2 "" -A 


. ^ 


(Vp) 


1/2 


(4.15) 


with the earth's magnetic field. 
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For infinite conductivity along magnetic lines, we have: 

E„ = 0 (4,16) 

and Ohm's law reduces to; 

Ei + vxBsO (4.17) 

Taking now the divergence of this equation and combining with (4.15), it 
is easy to arrive at: 

which relates the space charge, or potential difference, across field lines 
with the parallel current associated with the waves. (In our problem, space 
charge corresponds to the electi'omotive force applied by the tether between 
different field lines.) 

On the other hand, by combining Maxwell's equations: 

o c - 9- 

^ " 3t 

^ X B = Uq 

we arrive at: 

vv . E - v^E = Uq 1= (4.19) 

projecting along the magnetic field (z) direction and taking into account: 

H ^0 "W ’ 

Combining (4,20) with (4.18), we can write, for example; 

(v^ • Ep = Vy^ (V| . Ej_) (4.21) 

which tells us that the transverse space charge propagates at the Alfven speed, 
within the MDH framework and parallel to magnetic field lines. We will not 
apply these formulas to the tether-balloon system. Since the tether is itself 
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covered by an Insulator, the outside current flows only along the flux tubes 
cut through by the balloon and the Shuttle (*). To have an Alfvenic disturb- 
ance, we must refer to large balloon dimensions as indicated by (4,12), Then, 
supposing that the current is uniformly distributed across the balloon flux 
tube, we can write: 



(4.22) 


and, from (4.13) we obtain: 




1 

>^oVA 


Ej- 


(4.23) 


The perpendicular electric field, taking ohmic losses in the tether into 
account, is given by: 



%^Aw 

L 


(4.24) 


where is the total Lorentz field and (.'or motion perpendicular to B): 


Ej^o ~ VqB = 0.23 volt/m 


(4.25) 


is the tether's resistance and L is its length. 



lo 

1 ^ ir ^B 

2 f PqV^ 


Thus, we can write: 


(4.26) 


where: 


I. = rR 


1 


-io 


•0 2 B 

is the total parallel current in Alfven waves, which one would have for a 
perfectly conducting tether. 


(4.27) 


(*) Actually, the Alfvenic current flows at an angle 0^ with respect to field 
lines given by = arctan . In this case it is 0.55° and is not 
important for the following estimates. 
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Figure 23 reproduces the Mfvenic current as a function of balloon 
radius for a tether length L * 100 km and the two values of resistivity p » 
0,03jjS?m and p * 0,15ynm. We have taken v^ * 800 km/sec as an average value for 
the Alfv&n speed between 100 and 300 km of altitude, Note that the curves are 
really valid only for balloon radii Rg ‘i.> 20 sina (meters), according to (4,12), 
as the current calculated 1^^^ refers to waves In the hydromagnetic range of 
frequencies. We see from Figure 23 that the resistive limits to the current 
(1j^ s 6.21 amps and 1^ « 1,22 amps for the two respective resistivities) are 
not reached even at quite high values of the balloon radius. 

4 , 3 Comparison of Alfvenic Current with Current Due to Charged Particle 
Collection 

The current calculated in the previous section, for the Alfven 
disturbance associated with TSS motion, should now be compared with the current 
r in the tether, which results from particle collection from the end electrodes 

w 

(see Section 2.1), 

Figure 24 gives the ratio as a function of balloon radius, refer- 

ring to a passive system, configuration A, and two values of resistivity 
(obtained from the results of Figures 5 and 23), Although the comparison has 
validity only for large balloon radii, the Alfvenic current is always greater 
than the collection current. 

Assuming that the current in the tether system has to be carried along 
magnetic flux tubes through Alfven waves, the results of Figure 24 indicate 
that, for configuration A,, particle collection at the electrodes determines 
the value of the current in the system, In other words, the current in the 
Alfvenic wings (or, approximately, in the flux tubes cut through by the 
balloon and the Shuttle) is limited to the value 
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In configuration B, and with an electron gun at the Shuttle so as to 
keep the Shuttle potential low, we see from Figure 9 that substantially 
higher values of current are reached at relatively small balloon radii. 
We cannot compare such radii with corresponding values of because, In 
fact;, we cannot speak of an Alfvenic current. There remains, however, the 
possibility that this case might show a limitation to the current (other 
than the resistive limitation), which Is related to the capability of the 
electromagnetic disturbance produced by TSS to carry current away along 
magnetic lines of force. 


4.4 Physical Picture of the Alfven Wave System Associated with TSS 

When the Alfvenic disturbance in the large balloon's flux tube reaches 
the E layer of the ionosphere, it encounters a change in transverse conduc- 
tivity. The corresponding parallel current system then closes transversally 
to the magnetic field through Pedersen and Hall conductivities. A result 
of the dense ionospheric layer is reaction to the electric field of the 
wave is then a possible parallel reflection of the wave electric field. 

We can write: 


^up , „ jdown 


down 


(4.28) 


-up 


where the transverse electric field of the down-going wave, Ej'" is 
the electric field of the reflected wave, and the reflection coefficient R 
is given (Mallinckrodt and Carlson 1978) by: 

1 - X 


R = 


1 + X 


with; 




'Aw 


TT 0 A p 


(4.29) 


(4.30) 
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being the ratio between the Pedersen integrated conductivity i]p and the 
Alfven wave conductivity 

If X » 1, the electric field (and therefore the B plasma velocity) 

is zero at the boundary. The E layer of the ionosphere acts, in this case as 

a metallic boundary with frozen magnetic lines, On the other hand, if X « 1, 

we obtain at the boundary twice the amplitude of the incoming wave. In the 

latter case there is negligible conductivity of the ionospheric E layer, so 
that in the limit of a perfectly insulating E layer, the magnetic lines have 
no further identity in this region. 

To obtain numerical values for the reflection coefficient, we take as 


typical values for z (Hanson 1965): 

r 


p 


10 mhos 


at day time 


(4,31) 


3 X 10"^ mhos at night time 

Taking an average E layer altitude of h = 100 km, we can further use the 
following values of electron- density : 

'2 X 10^ cm"^ during the day 


e 


2 X 10^ cm-^ during the night 


(4,32) 


to determine the Alfven speed at the E layer. The final, typical results for 
the reflection coefficient R are: 

-0.92 during the day 




+0,17 during the night 


(4.33) 


Thus, whereas there is almost complete reflection from the day-time ionosphere, 
reflection is quite weak from the night-time ionosphere, and most of the 
Alfven disturbance is transmitted further downward. 

The next point to consider is the propagation time of Alfven waves from 
the tether system to the E layer levels. This time varies, of course, from 
one end to the other of the tether causing reflection (or partial reflection) 
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of different parts of the wave front at different times. 

Taking a typical distance of 100 km and an average Alfven speed of 
800 km/s, we arrive at a round trip time of the Alfven wave; 

T 'V 0.25 sec 

in vffiich time the tether balloon system has moved by a distac.ce: 

Ay 'V 2 km 

Clearly, therefore, when there are reflected day-time waves they will no 
longer find the tether-balloon system on their way back from the E layer. 

Thus, an influence of the ionospheric E layer on the current in the tether 
system, which was indicated in the case of the moon of Jupiter lo (Goldreich 
and Lynden-Bell 1969), is not possible. 

Correspondingly, the picture of a so-called "dc current circuit" moving 
with the tether— i^e. , a freezing with the tether's motion of the intercepted 
flux tubes (which would then slip with respect to the ionospheric base)-- is 
not appropriate. 

What happens in daytime is that the Alfven waves radiated by the tether 
are partially reflected from the ionospheric E layer. They then travel back 
to the conjugate ionosphere, where they are again partially reflected, and 
the process repeats for a number of times. During its motion, the tether 
system generates a system of waves reflected back and forth betv/een conjugate 
zones of the E layer, all along its orbit. 

It is interesting to ask how many reflections are possible before the 
amplitude of a given wave decreases significantly. When reflections occur 
with the same reflection coefficient r , the number of successive reflections 
necessary to achieve a ^ reduction of the wave amplitude would be: 

^ ^ ^ ^A 

A 
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Thus, N 'V 12 for reflections between daytime ionospheres. In the case of 
nighttime ionospheres (r 'v- 0.17), more than 80?o of the Alfvenic disturbance 
is transmitted downward. 


4.5 Power into Guided Alfveri Waves 

Alfven waves are guided along the magnetic tubes of force, i . e . , their 
energy is propagated along magnetic field lines. As a consequence, the power 
in Alfven waves can be obtained by multiplying the wave energy density W for 
the volume filled up by wave energy in a second, This volume, in turn, will oe 
given by the cross-section of the flux tubes intercepted by the system (A) 
multiplied by the wave group velocity Vq. 

Thus: 

= W X 2A X Vg (4.35) 


where the factor 2 indicates wave propagation in two opposite directions with 
respect to the tether system (down to conjugate regions of the low ionosphere) 
For Alfven waves, it is easy to show that: 

_ b2 


W 


Aw 2Ur 


(4.36) 


being the magnetic field of the perturbation. That is, the energy of the 
fluctuations is entirely magnetic. Thus, again using Vg = v^^, we obtain for 
the power P^^^ associated with Alfven waves: 




A v^ 


(4.37) 


As it refers to Alfven waves, this formula is valid only for sufficiently 
large dlriienslons of the conductor in the direction of motion [see (4.12)]. 
For the TSS system, it will, therefore, apply to power contained in the bal- 
loon cross-section, where the balloon radius satisfies: 
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Rg > (20 sina) meters 


(4,38) 


In terms of the flux tubes cut through by the conducting part of the Shuttle 
and the tether itself, the formula is not strictly applicable. This is because 
the corresponding electromagnetic disturbances contain components of frequencies 
above the hydromagnetic range, which neither have an energy density given by 
(4.36) nor propagate parallel to _B with group velocity v^^. We will, therefore, 
use (4.37) to calculate the Alfvenic power specifically in the balloon flux 
tube. 

The perturbation magnetic field in (4.37) is calculated in terms of the 
transverse electric field across the balloon cross-section; 


B-viE, 

and, taking into account resistive losses in the tether: 






(4.39) 


(4.40) 


wi th : 


!l.0 ' Vo 


(4.41) 


Thus, we obtain for the Alfvenic power: 

n P,2 

Pnu TtRr 

AW B 


(4.42) 


In (4.42), referring to configuration A, we will have to use the current; 

■i = i'c(Rb) 

due to charged particle collection (see Section 4.3). 

An input impedance of the balloon flux tube, with respect to guided 
Alfven waves, can now be easily calculated from: 


Av" 


"AW p 


(4.43) 


AW 
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where: 


AV 2RgE^ 


is the potential difference across the flux tube. From (4.42) and (4.43) we 
obtain: 


Z 



^o^'a 


(4.44) 


A typical value of this impedance, at the altitudes of interest to us, is: 

-V 1 . 3 ohms 


Figure 25 shows our results for the power '"aw* as a function of balloon radius. 
We are referring to configuration A, a passive system, and the two values of 
resistivity p = 0.03yf2m and p = O.lSynm. For example, for Rg = 40m, we obtain: 


^ 200 watts 
-V 60 watts 


for p = 0.03ynm 
for p = O.lBynm 


The difference between the two curves, and, in particular, the knee in 
curve b) that refers to the case of higher resistivity, are understandable in 
terms of the concepts that we have already introduced and of the formula^ (4.42) 
that was used to compute P^^^. Curve b) reaches a plateau between Rg = 50 m 
and Rg = 90 m. This is due to the fact (see Figure 5) that, for the same 
radius, the current, and hence the ohmic losses, that appear in equation (4.42), 
do not vary very much. At about Rg = 90 m, the balloon potential drops to a 
small value (see Figure 4) and, correspondingly, the Shuttle potential 
increases. Hence the current start increasing again (see Figure 5), and this 
determines an increase in ohmic losses. This more than balances the increase 
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of with Rg and explains qualitatively the decrease in for Rg- s 90 n. 

For Rg > 100 m, increases again. This last trend is due to the fact, 
already mentioned in connection with Figures 17 and 18, that for p = O.lSynm 
and Rg > 100 m, the current, and hence the ohmic losses, reach an asymptotic 
value. 

The electric field, then (see equation 4.40), does not vary any more 

2 

and we see from the same equation that •"aw starts increasing as Rg. This 
behavior does not appear in curve a) of Figure 25. It would appear at much 
higher values of Rg, not shown in the Figure, owing to the smaller value of 
the wire resistivity, and hence to the reduced relevance of the ohmic losses. 

The reader should take note that what we have calculated here is strictly the 
power contained in the flux tube cut through by the large balloon. The tether 
enters the calculation only through the ohmic reduction of the transverse 
electric field across the tube [see (4.40)], due to its internal resistance. 
This power has, therefore, nothing to do with the total power radiated by TSS 
as a driven antenna (if its current were pulsated) at different frequencies. 

It is in this last power that the radiating properties of the tether as a long 
dipole (and, therefore, its length) will enter in an essential way. 

4 . 6 Preliminary Considerations on Wave Detectability on the Ground 

4.6,1 General 

A source of electromagnetic radiation in the ULF/ELF band, embedded in 
Earth's ionosphere (Sypi cal height 200 km), illuminates Earth's surface through 
several possible propagation mechanisms. One of them could be properly called 
the "artificial micropulsation" mechanism. This mechanism consists of Alfven 
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waves that are guided by the lines of force of Earth's magnetic field, from 
satellite height to the bottom of the ionosphere, spread from there downward 
semi -spherically and directed to Earth's surface as horizontally polarized 
downward e.m. waves similar in many respects to a natural micropulsation. 

An estimate of the expected field intensity for this case can be formulated 
on the basis of the geometric optics treatment worked out by Kelly et al . 

(1974, 1976). This is based on the computation of the "wave-spreading factor," 
as illustrated in Figure 26. 

A second mechanism consists in the excitation frequencies at ULF/ELF 
of the Earth-Ionosphere wave guide from a satellite-borne antenna. This 
mechanism was analyzed by Einaudi and Wait (1971a, 1971b) for the case of a 
flat Earth and of a vertical geomagnetic field orientation, and by Pappert 
(1973) for spherical Earth and arbitrary magnetic field orientation. These 
authors have all considered an infinitely small dipole as the satellite-borne 
radiator. With this second mechanism, the satellite-borne antenna excites 
in principle, long-range, quasi-TM and quasi-TE modes, which propagate 
horizontally in the Earth- Ionosphere wave guide, away from the sub-satellite 
point. Some of these modes are characterized by very low attenuation. How- 
ever, the excitation factors involved are exceedingly small, especially for 
the case of the vertically polarized satellite-borne elementary dipole," 
which is of direct interest here. 

4.6.2 Estimates of E.M. Field Intensitives at Earth's Surface 

By following the geometric optics approach of Kelly et al . (1974), we 
compute first the radiation intensity U (power per unit solid angle) of the 
orbiting antenna. In the case of an isotropic radiator, with output power 
P, we have: 
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Figure 25. Power In Alfven wave from a single balloon 
versus balloon radius! Configuration A, 

(a) 0.03 


(b) ^ -OjlS^uiin, 



Figure 26. Spreading factor quantities diTL and. 

dAg. The multihop ares element dAg^ ^ 
is also shown. 

( from Kelly et.al., 1974) 


In our case, which pertains to the excitation, by the orbiting antenna, of 
Alfven waves, we have: 



because all the power is channeled along the line of force of the geoinagneti : 
field passing through the source, with half of the power propagating upward 
and half downward. The power density P, at Earth's surfaoe (vectors 

E and H are horizontal) is then obtained by multiplying the radiation intensity 
U by the spreading factor ^ 1 ^ (in”^), computed with the geometric optics 
approach: 

p _ ,1 1 dQ 

^d " ^ HBsI ’Sq 

Examples of the computation performed by Kelly et al . (1974) are given in 
Figures 27 and 28. These two figures define the symbols used in the expression 
of the spreading factor and give its dependence upon latitude and longitude 
for the frequency of 3kHz, for two altitudes of the orbiting antenna, and for 
three models of the ionosphere. Model Al is characterized by a linear increase 
of electron density from 70 km to 300 km, with TO® el/cc at 70 km and 2 x 10® 
el/cc at 300 km. Model A2 is characterized by a realistic profile of the 
bottom-side ionosphere, with 10® el/cc at 84 km, 10** el/cc at 200 km, and 
2 X 10® el/cc at 300 km. Model A3 is characterized by a latitudinal dependence 
and by a smoothly varying electron density with height. All these models are 
described in the Appendix to Kelly et al . (1974), 

In order to compute, without performing ray tracing, the spreading 
factor at frequencies in the ELF/ULF band, other than the VLF ones chosen 
by Kelly et al . (1974) in their examples, the following approximate equality 
can be used: 
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spreading factor - ^ 

where \i is the Index of refraction of the wave at satellite height and h is 
the height of the bottom of the ionosphere above ground. Table HI taken from 
et al . (1974) compares rigorous and approximate calculation of the 
spreading factor, for the frequency of 3kHz, for several satellite heights, 
for various latitudes and for three ionospheric '">uels. 



LATITUOF. (deoreets). 


Figure 27. 


Variation of spreading factor with 
receiver latitude for two fixed 
transmitter locations. The X's and 
O's are results for lower ionospheric 
Model A-2. The linos are results for 
Model A-1 (from Kelly et al . , 1974). 


SPREADING FACTOR 



Figure 28. 


Variation of spreading factor with 
receiver longitude for two fixed 
transmitter locations. Model A-1 
results are shown (from Kelly 
et al .. 1974). 
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Table III 


Comparison of the 3-kHz spreading factor obtained by geometric optics calcu- 
lation (-J— A), with the approximate value (-ir). 
cosa dME jrh^- 


Transmitter Coordinates 
Height (km), Latitude 

Lower Ionosphere 
Model 

Approximate Value* 
n*’h^' 

Geometries 
Optics Ex- 
pression 

-JL-. (km“^l 

cosa dAg ^ 

818 

43.6»N 

A-1 

6.8 X 10“^ 

4.3 X 10“^ 

818 

43,6'»N 

A-2 

4.7 X 10-? 

2.8 X 10"^ 

818 

43,6'’N 

A-3 

4.6 X 10-^ 

X 

o 

1 

2152 

38.5<’N 

A-1 

8.95 X 10-7 

2.0 X 10-7 

2152 

38.5°N 

A-2 

6.21 X 10~7 

1.6 X 10-7 


* For comparisons to ray-tracing prediction using Models A-1 , A-2 
and A-3, heights h of 70, 84, and 70 km, respectively, In the 
expre/islon . (Ionospheric Models A1 , A2, A3, as described 
in the Appendix to Kelly et al . « 1974). 





Table in shows that, at least for the computation conditions adopted, 

the approximate value of the spreading factor is accurate enough the for pur- 

1 

pose of our study, By assuming that the accuracy of the expression for 

u ’h 

the spreading factor holds true at ULF/ELF frequencies, we have the following 
value of power density at Earth's surface and at the frequencies of interest: 


P 


d 




where P is assumed to be 100 watts. 

M = with c « velocity of light in free space « 3 x 10^ m/sec 

V = velocity of Alfven waves = 8 x 10^ m/sec 

h = height of ionospheric bottom above Earth's surface » 6 x 10'^ m 


Therefore: 

p =, 100 1 

d ' 2 ' (375)M6 x10")^ 


10“ watts/m^ 


Hg = 1.62 X 10-® ^ (-156 db with respect to ^). 

Eg == 6.14 X 10-^ I (-104 db wrt 

(both Hg and Eg are horizontally polarized). 

The noise levels that the signal has to overcome are: 

- internal noise of a cryogenic magnetometer, = 7.9 x 10"^ At/m 

(-162 db//Hz with respect to this is -172 db wrt 
when using 0.1 Hz bandwidth) 

- external micropulsation noise (discontinuous occurrence) 

H[^ = 7.9 X 10"^ ^ (-142 db//Hz wrt i^), for horizontally 

*1 f' 

polarized This is -152 db wrt — rjen using 
0,1 Hz bandwidth). 


and 


,y 
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Therefore, the signal ■‘to-noise ratio Is approximately +16 db (In 0,1 
Hz band-width) In the time Intervals during which micropolarization activity 
(for Instance, "pearls," with carriers at a few Hertz) subsides, and the only 
noise to contend with It the cryogenic magnetometer's Internal noise. On the 
contrary, the signal-to-nolse ratio (SNR) Is approximately -4 db (still In 
0.1 Hz bandwidth) when micropul satlon activity occurs. We think that this 1s 
a.lequate enough for an Initial experiment, considering that mi cropul satlon ac- 
tivity 1n the frequency band of interests (ELF/ULF) 1s a discontinuous phenomena. 
Besides, in such a case as natural "pearls," the trains are made of bursts 
(for Instance, of 1 c/s carrier) alternating with pauses, which makes It 
possible to detect the satellite-emitted signals even with an SNR as low as 
-4 db. Signal Integration with longer integration times would further enhance 
detectability. 

Let's consider now the case of the excitation of quasi -TM and quasi -IF 
guided modes in the Earth-ionosphere cavity, due to an infinitesimally small 
elementary dipole in orbit. Limiting ourselves to an example with a carrier 
frequency of 75Hz (the frequency of the numerical example worked out by 
Pappert, 1973), we have the same field intensities as Pappert if we choose 
an orbiting antenna 318 km long, with a current of 10 amp (electric moment = 

3.18 X 10® Am). The signal intensities to be expected at Earth's surface can 
be read in this case along the y-axis scale, which is to the left in Figure 
29 and along Earth's surface of several niegameters from the vertical that 
contains the orbiting radiator. However, the tether antenna presently under 
consideration by the designers of the TSS facility has a length of 100 km. 

If the current acquires the realistic value of 1.79A (electric moment - 
1.79 X 10® Am), the signal intensities are expected to be smaller. Their 
values can be found by changing, in Figure 29 and Figure 30, the scale of 
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NOTE 

Dip is at 0° at the 
magnetic equator 


Figure 29 


Sir,nal levels versus distance, Signal levels for eleetric- 
dipole sources referred to a current moment of 3.18 X 10* 
amp-m. Signal levels for magnetic dipole sources referred to a 
current loop of 2.02 X lO'* arnp-m’. The nzlmulh is 90“, the 

dip is 75°; and the frequency is 75 Hz. Legend: 

ground-based electric dipole, end fire, a «=■ 10"* mho m“‘; 
• c • • • 4 vertical electric dipole, 500 km; hori- 
zontal electric dipole, broadside and end fire, 500 km; 

horizontal magnetic dipi broadside and end fire, 500 km. 

(from Pappert, 1973). 
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+15 

+10 

+ -5 ' ’ 

^ dB v;ith respect 

-5 to 1 ^iV/m 

-10 

-15 

-20 

-25 


NOTE - Dip is 
.0 

0 at magnetic 
equator 


Figure 30 

Signal levels versus distance. Signal 
levels for electric-dipole sources 5 
referred to a current moment of 3.18 10 
amp.m. Signal levels for magnetic dipole 
sources referred to a current loop of 
2,02 10 ^amp.m^. The azimuth is 90°, the 
dip is 15 , and \,ne frequency is 75 Hz. 

Legend : — ground-based electric 

dipole end fire, = 10 ”^ mhos/m; 

• ••»<?•# vertical electric dipole, 

500 Kmj — horizontal electric 

difsole, .end fire, 500 Kmjooooooo#. 
horizontal electric dipole, broadside, 500 Km ; 
y \ y X X horizontal magnetic dipole, end 

fire, 500 Km; horizontal magnetic 

dipole, broadside, 500 Km. 

(from Pappert, 1973) 
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the y axis from the one on the loft side (Pappert, 1973) to the one on the 
right side. This decrease of the electric moment causes, in fact, a 25 db 
decrease in signal intensity. Figure 29 indicates that for a dip of 75°, 
and at a distance of 500 km (along the surface of Earth) from the vertical 
that contains the satellite (elementary dipole, vertically oriented, at a 
height of 500 km above Earth's surface), the signal intensity is approximately 
-7.5 db wrt (= 0.42 ~). Under the same conditions, except with the dip = 
15°, Figure 30 gives a signal intensity of +7.5 db wrt (= 2.38 ^). If 
we observe with a cryogenic magnetometer the horizontally polarized magnetic 
field of the wave (integration time adopted = 10 seconds, as before), we have 
a SNR that varies between +8 db and -7db for a dip inclination's change from 
15° to 75°. This assumes that the only noise we have to contend with is the 
cryogenic magnetometer's internal noise. 

It would appear that the illumination of a "hot spot" (Kelly et al . . 
1974) on Earth's surface, by simulating an artifical micropulsation, is a 
more favorable mechanism (at least at close range from the vertical that 
contains the orbiting radiator) than the excitation of quasi -TM or quasi -TE, 
long-range, guided modes in the Earth-ionosphere cavity. 

This conclusion holds for the case in which the model adopted for the 
orbiting radiator is an infinitesimally small elementary dipole, (but it is 
characterized by the total moment of the long wire), and is vertically 
oriented. We recommend possible follow-on efforts that adopt the more 
realistic, albeit more complex, model of an orbiting antenna with its full 
length (100 km or larger). 
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5. 


Conclusions 


The report has investigated the potential capability of the electro- 
dynamic tether to generate electromagnetic waves in the ionospheric plasma. 

We presented several results concerning current-voltage characteristics of the 
tether. It has been pointed out that, in a purely passive configuration-- 
either with the tether deployed downward or upward, relatively large potentials 
may develop at the Shuttle. To avoid this, one should consider using charged 
particle guns. The calculations have included the effect of secondary electron 
emission from the balloon surface, which can be significant for configuration 
A when the balloon is at high potentials with respect to the plasma. The 
final scope of these calculations indicates the possibility of navir;g large 
currents in the tether and of computing, correspondingly, the power dissipated 
in ohmic losses, the power going possibly into accelerating secondary electrons, 
and the remaining electrodynamic power that 's, in principle, available for 
wave generation. 

An important conclusion is that significant powers are available for 
wave generation, even when effects of secondary electrons are taken into account 
(and in the framework of an overestiniation of that effect). For configuration 
A (tether deployed downwards), we find powers of the order of some kilowatts 
even for balloon radii as small as '^'5 meters, and greater for larger balloon 
radii. For configuration B (tether deployed upwards), if one uses an electron 
gun at the Shuttle in such a way as to keep it at plasma potential, one obtains 
both large currents and correspondingly large powers for wave generation at 
very small radii. For example, for Rg = 2 meters, the current is i 'v 4 amps 
and the power is above 10^ watts. 
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In the case of a particular electrodynamic phenomenon and a program 
for observing it, (either from a sub-satellite in the ionosphere or from the 
ground,) the above calculations of current-voltage characteristics constitute 
the basis on which to decide what parameters of the TSS system are necessary 
to make that phenomenon detectable . 

A consideration of the possible effects generated by accelerated 
secondary electrons led to the following conclusions: 

1) In terms of the velocity distributions of secondaries outside the 
sheath region surrounding the balloon, excitation of electromagnetic instab- 
ilities (requiring an excess momentum perpendicular to the magnetic field) 
is not likely to be Important from an energetic point of view. 

2) A phenomenon that may possibly be important and that merits further 
investigation is the possible excitation of UV radiation by high-energy 
electrons moving down the balloon flux tube and, in the absence of other inter- 
interactions, reaching the Earth's atmosphere. This conclusion may be drawn 
from a comparison of calculated power densities in energetic secondary elec- 
trons (3 X 10^ ergs car^sec-ifor a balloon radius Rg = flm) with corresponding 
numbers referring to the actual measurements of Lyman-a/UV excitation at the 
foot of the lo's flux tube. 

The next topic considered is what we have called "spontaneous" wave 
generation, i^. , the electromagnetic disturbance that is naturally associated 
with the motion of TSS across Earth's magnetic field. In particular we have 
focussed our attention on the possible spontaneous generation of low-frequency 
Alfven waves. This requires rather large balloon dimensions (Rg _> 20 meters, 
for motion perpendicular to B) and power levels that one gets for Alfven waves 
guided in the flux tube cut through by the balloon are significant: > 100 

watts. We have also determined the feasibility of detecting this Alfvenic 
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disturbance on the ground, through the use of cryogenic niagnetometers . 

As we have discussed, to have significant power in spontaneously guided 
Alfven waves, one needs large balloon dimensions. On the other hand, by 
modulating the natural current in the tether (through modulation of the 
electron gun at the Orbiter in configuration B), one would use TSS as a 
long antenna and then take advantage of the properties of the long tether 
itself as a radiator. This possibility deserves further study. 


« 
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